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The propagation and transmission of Bessel beams through nano-layered struc- 
tures has been discussed recently. Within this framework we recognize the for- 
mation of unguided diffraction-free waves with the spot size approaching and 
occasionally surpassing the limit of a wavelength when a Bessel beam of any 
order n is launched onto a thin material slab with grazing incidence. Based 
on the plane- wave representation of cylindrical waves, a simple model is intro- 
duced providing an exact prescription of the transverse pattern of this type of 
diffraction-suppressed localized waves. Potential applications in surface science 
are put forward for consideration. 
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1. Introduction 

Nondiffracting Bessel beams are associated with solu- 
tions of the Helmholtz wave equation in free space, which 
transverse amplitude distribution is written in cylindri- 
cal coordinates as the product of a linear phase function 
exp (inO) depending on the azimuthal coordinate and a 
radially- varying Bessel function of the first kind and or- 
der n 1 . The maximum value of the zero-order Bessel 
function is attained at the origin; however, higher-order 
Bessel functions have a phase singularity at r = 0, reach- 
ing its highest intensity nearby 2 . Its significance relies on 
the possibility of engineering laser beams reaching high 
intensities in a reduced area without diffraction-induced 
blurring. In particular, the envelope of the Bessel func- 
tion falls off as 1/ y/r, which leads to an electromagnetic 
field having an intensity distribution that is not square 
integrable 3 . As a consequence finite-energy practical re- 
alizations of Bessel beams show axial bounds between 
which diffraction-free propagation is valid^*2*Z*£. 

The propagation and transmission of Bessel beams 
through nano-layered structures has been discussed 
recently ^ 10 ! 11 ! 12 in the context of superluminality in op- 
tical wave packet o 13 i 14 i 15 i 16 i 17 . However, an important 
attribute such as diffraction-free propagation is main- 
tained only if the optical axis of the beam is parallel 
to the medium interface ^ 18 ! 19 . Guided waves in planar 
linear 20 and nonlinear— waveguides are of this kind of 
modal fields having exponentially decaying tails far be- 
yond the core layer. In general this condition relies on 
grazing incidence upon the interfaces, which may be used 
to study surfaces 2 ^ by increasing wave penetration. 

In this paper we conveniently introduce a simple model 
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from which high-order Bessel beams are interpreted as 
a coherent superposition of propagating and counter- 
propagating waves along a direction (normal to the layer 
interfaces). In fact this sort of cylindrical wave local- 
ization is essentially understood as the generation of a 
line focus around the z axis with 2tt illumination. We 
recognize the formation of diffraction-free waves when a 
Bessel beam of any order n is launched onto a thin ma- 
terial slab. The resultant optical beam has in general a 
twin-peak pattern, which spots size even surpasses the 
limit of a wavelength. From a geometrical approach the 
thin film gives rise to a pair of line images associated with 
waves impinging onto either the left or the right face of 
the layer. 

2. Vector Bessel beams 

The electric and magnetic fields of a monochro- 
matic diffraction-free wave travelling in an homogeneous 
medium are written using the modal function ansatz 

E(x, y, z, t) = e(x, y) exp (i/3z — icut) , (la) 
H(x, z, t) = h(x, y) exp (i/3z — iuot) , (lb) 

where uo is the time-domain frequency of the radiation 
and (3 > is the propagation constant along the z axis. 
We consider a non-absorbing medium with positive di- 
electric constant (e > 0), zero magnetic susceptibility, 
and free of charges. Let us assume a polarized wave 
such that the ^-component of the electric field vanishes, 
E y = 0; complementary polarized waves with H y = 
may be obtained straightforwardly by substituting fio <-> 
eoe, E — > H, and H — > — E. Faraday's law leads to the 
estimation of the magnetic field in terms of the electric 
field, H = (—i/ujfio) V x E. Furthermore we obtain the 
equation E z = i(3~ 1 (dE x /dx) from the Coulomb's law. 
Therefore we may use the two-dimensional scalar func- 
tion U (x, y) = e x in order to describe the wave mode un- 
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Fig. 1. Normalized amplitude (upper row) and complex 
argument (lower row) of the wave function U n (x,y) for 
ft = 9.07 /im" 1 and indices (a-b) n = 0, (c-d) n = 1, and 
(e-f) n = -2. 



ambiguously. Let us point out that alternate routes for 
the description of electromagnetic diffraction-free beams 
may be found elsewher o 23 i 24 i 25 . 

The function U(x,y) satisfies the differential wave 
equation 



d 2 U d 2 U 



dx 2 dy 



(2) 



where the wavenumber k = co^e/c and c is the speed of 
light in vacuum. Particular solutions of the wave equa- 
tion ([2]) that are periodic in the azimuthal coordinate 
U n = / n (r) exp (inO) being r the radial coordinate and 
n an integer, leads the radial function f n to satisfy the 
Bessel's equation 1 . Assuming / n (0) < oo and (3 < k we 
may write f n (r) = J n (ftr), where J n is a Bessel func- 
tion of the first kind and ft 2 = k 2 — (3 2 . In Fig. [1] we 
plot the 2D wave field U n of different orders n at a wave- 
length A = 0.6 jim (k = 2tt/X) and propagation constant 
P = k/2 (K, = y/3/3). 

3. Radiation modes with on-axis localization in- 
side a thin layer 

In this paper we consider the monochromatic non- 
diffracting beam is propagated in the layered medium 
of Fig.Efa). The y axis is set such that it is perpendicu- 
lar to the surfaces separating the centered slab (of width 
L) with dielectric constant 



e(l + S), 



(3) 



and the cladding semi-infinite layers of dielectric con- 
stant e. Therefore, E y = (H y = 0) indicates that the 
diffraction- free beam is s-polarized (p-polarized). The 
existence of the Bessel beam is assumed in the region 
| y | > L/2. In the limit L — > it obviously yields a 
radially-symmetric perfect Bessel beam. Our purpose 
is the analysis of the Bessel-beam transformation in the 
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Fig. 2. (a) Schematic geometry of the planar-layer-based 
medium. In (b) we represent the wave fields in regions 
established by interfaces at y = ±L/2. Perfect matching 
e r = e leads to wave fields U ± = F ± = R ± . 



presence of the thin film. For the sake of clarity in our 
formulation, let us first consider the limiting case setting 
either L = or S = (perfect match, e c = e). 

A. Propagating and counter-propagating modal decom- 
position 

From a practical point of view it is necessary to describe 
the wave field by means of propagating and counter- 
propagating fields, U n = + U~ . The propagating 
field [/+ is assumed to be excited at y — > — oo travers- 
ing the semi-infinite layer in the positive direction of the 
y axis, as depicted in Fig. E^b). On the other hand, 
U~ emerges from y — > +oo to go through the right-side 
cladding layer directed in the negative direction. One 
can see that an approach describing U n in rectangular 
coordinates is preferable. 

To obtain a convenient representation of the fields 
and U~, let us recall the Sommerfeld's integral represen- 
tation of the Bessel functions 1 , 



Jn(p) 



exp (zpcosc 



(4) 



Consequently, a Bessel beam can be cast into homoge- 
neous plane waves whose wave vector directions form a 
circular cone around the z axis 26 , 



2tt 



exp [ina + in (x cos a + y sin a)] da. 

(5) 



Denoting 7 the angle measuring the aperture of the cone 
we obtain n = ksiwy and (3 = k cos ^y. In Fig. [T] 7 = 
60 deg. When < a < 7r, the projection of the wave 
vector onto the y-axis ft sin a is positive, representing a 
wave field whose phase advances toward positive values of 
y\ this shall be written as . Otherwise ft sin a < for 
— 7r < a < leading to a counter-propagating field, U~ . 
Setting k x = ft cos a in Eq. ([5]) we obtain a ID Fourier 
representation of the field components 

1 f K 

Un(x,y) = — / U n t (k x )exp(±ik y y)exp(ik x x)dk x , 



2tt 



(6) 
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where the spatial spectrum 

rhy rv 

and k y = n 2 — k 2 . 

We point out that the cylindrical wave Uq given in 
Eq. ([5]) for n = may be essentially understood as an 
aberration-free 2D focused field with focus at the origin, 
x = y = 0, using 2tt illumination. Setting net = k$> 
allows us to identify the aberration function^ <I> that 
would depend linearly upon the Bessel index n and the 
azimuthal coordinate a. 

Let us conclude our analysis on /7+ by giving some 
functional expressions for its fast computation. For non- 
absorbing dielectric media (e > 0) we obtain the relation 

U~ = (-!) n (^n )*• This means that 

U;{x,y) = {-l) n [U+(-x,y)]*, (8) 

so that propagating and counter-propagating compo- 
nents of the Bessel wave functions possess mirror sym- 
metry with respect to the plane x = 0. Furthermore, 
using the Jacobi-Anger expansion^ in Eq. ([5]) we may 
obtain the analytical expression 

CO 

U+{x,y) = - smc[(m-n)/2]U m (x,y), (9) 

171 — — CO 

where the function sinc(u) = sin(7ra) / (ttu) and sinc(0) = 
1. In Fig. [3] we represent U+(x,y) associated with the 
Bessel beams shown in Fig. [TJ Radial symmetry char- 
acteristic of Bessel patterns is not conserved here, ex- 
hibiting an elongated peak along the y axis. Moreover, 
this focus is shifted toward positive (negative) values of 
the x coordinate for n > (n < 0). Finally, the on- 
axis vortex^ for the order n = — 2 is maintained (not for 
n = 1); however in Fig.(3](f) we observe that the topolog- 
ical charge has changed. Additional phase singularities 
come out at points all along the x axis even for n = 0. 

B. Transverse wave fields 

Let us now evaluate the field transformation of the Bessel 
beam U n in the presence of the core layer. We assume a 
wave field J7+ [see Eqs. (j6]) and Q] propagating in the 
semi-infinite medium y < —L/2 and impinging at the 
core-cladding left interface, as seen in Fig. Efb). Also a 
counter-propagating field U~ evolving in the semi-infinite 
medium y > L/2 is incident at the right face of the core 
layer. In general we may represent the x-component of 
the electric field transmitted in \y\ < L/2 employing a 
factorization of the form given in Eqs. ([T]). Also sep- 
arating the propagating and counter-propagating parts 
of the field as F+ + F~ we may use the Fourier inte- 
gral representation shown in Eq. (j6]) substituting k y by 
k c ^ y — ^k 2 — /3 2 — k\ , where k c = ujyJT c /c. Explicitly 
we write 

1 f K 

F±(x,y) = — / F±(k x )exp(±ik Cjy y)exp(ik x x)dk x . 

(10) 



n = n = 1 ?7 = -2 
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Fig. 3. Wave function U^(x^y) associated with Bessel 
beams of Fig. [U Similarly, normalized amplitude is rep- 
resented in the upper row and complex argument is plot- 
ted in the lower row. 



Here k c , y = k' c y + ik" y is in general a complex number, 
where k' c > 0; also k" > leading to either homoge- 
neous or damped waves inside the slab. 

Additionally, the fields R~ (in y < —L/2) and 12+ 
(in y > L/2) flying away from the core layer should be 
included for a complete description, schematically shown 
in Fig. EJb). In core barriers, these fields account for 
reflected waves from and U~ , respectively. However, 
in the limit 5 — > we would have 12+ = related to 
transmitted waves through the core medium. In general, 

Rt(x,y) = — / Rt(k x )exp(±ik y y)exp(ik x x)dk x: 

(11) 

stand for a combination of reflected and transmitted 
waves as shown below. 

We point out that the condition /3 < k get rid of a 
diverging behavior of far off the core layer. These 
onedirectional wave fields on other hand lead to the for- 
mation of localized radiation modes. Such an optical feed 
is necessary to obtain a real- valued propagation constant. 
Otherwise it gives rise to leaky modest being (3 a com- 
plex number, which might be understood as disability to 
suppress the natural intensity decay induced by diffrac- 
tion. 

In the following, the spatial spectra F+ and 12+ are 
estimated. Boundary conditions for the s-polarized waves 
at the core-cladding interfaces lead to the continuity of 
the scalar field E x and its normal derivative [i^ x dE x /dy 
at y = ±L/2. Therefore the propagation constant /3 is 
conserved in the layered media. Also we obtain 
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where 
A = 

B = 
C = 
D = 



ex.p[i(k c ^ y — k y )L] 
{k c ,y ky) — (k>c,y ~ ky) exp(2i/c c?: y_Zv) 
2i { k c, y ~ k l) sin(fe c>3/ L) exp[i(k c ,y - k y )L] 
(kc,y H~ ky) {kc,y ky) exp(2zfc c?: yZ/) 

'Iky {k Cy y + ky) CXp [ft ( /C c ? y ~ ky)L/2] 



(13a) 

, (13b) 
(13c) 



(kc,y H~ ky) (k c ,y ky) exjp(2,ik c ^yL) 
2k y (k Cj y - ky) exp(ik c , y L) exp[i(k Cj y - k y )Lj2^ 
(kc,y H~ ky) (k c ^y ky) exjp(2ik c ^y L) 



= ky. How- 
k y \L > 1 for 



Excepting the phase factor exp (—ik y L), the terms A(k x ) 
and B(k x ) are simply the Airy's formulas corresponding 
to the transmittance and reflectance coefficients through 
the slab2£. Thus, the matrix elements A and B quantify 
the contribution of reflected and transmitted wavelets to 
the scattered fields R^. 

We confirm that = F^ = under perfect match- 
ing conditions; note that 5 = yields k CyV 
ever, serious deviations arise when \k c ^ y — 
which the phase (exponential) factor in the numerator 
of the coefficients ABCD given in Eqs. ([T3]) oscillates 
(falls off) rapidly. This inequality is approximately sat- 
isfied if \S\L > 2n/k 2 . Specifically for both material and 
beam sizes in the order of the wavelength, L ~ A and 
Kj ^ /c, a case that is analyzed in the next section, it 
yields \5\ > 1/tt so that the difference of refractive in- 
dices for the core and the cladding layers might not be 
large in excess. 

4. Numerical simulations 

Let us first evaluate numerically the amplitude distri- 
bution of the wave function e x corresponding to Bessel 
excitations as those of Fig. [1] in a centered dielectric slab 
of dielectric constant with relative increment 5 = 10 -2 . 
Therefore the source spectra are those of Eq. ([7j). 
We choose L = 100 jam in order to hold the inequality 
\S\L > (sin7)A/7r at A = 0.6 \im and 7 = 60 deg. In this 
case the slab width is much larger than the spot size of 
the Bessel source so that wave localization is expected at 
points far from the surfaces. 
The resultant wave function, 

U+(x,y)+R-(x,y) ]£y<-L/2; 
e x (x,y) = { F+(x,y)+F-(x,y) if \y\ < L/2; (14) 
U-(x,y)+B+{x,y) if y > L/2, 

is depicted in Fig. 2J Irrespective of the Bessel index n, 
the central high-field patterns shown in subfigures from 
the upper row consist of twin spots with large sidelobes. 
For n = the spots appear along the y axis but an addi- 
tional lateral displacement may be observed in the trans- 
verse x direction if \n\ > 0. Also a plethora of phase 
dislocations are found in all cases. 

The geometrical approach (A — > 0) bring us a sim- 
ple interpretation for the results of this simulation: We 
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Fig. 4. Transverse distribution of Bessel-excited wave 
field e x (x, y) inside a thin slab of width L = 100 jam made 
of a material with relative dielectric constant 5 = 10 -2 . 
Again \e x \ is plotted in the upper row and its complex 
argument is shown in the lower row. 
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Fig. 5. Ray tracing representation of the Bessel-like 
beams shown in Fig. [H Dual parallel axes at y = 
±0.25 fim (and x = 0) set the confinement region of the 
nondiffracting beam. Spherical aberration drawing the 
envelope of the refracted rays gives rise to a bilateral- 
symmetric caustic. 



may think of the incident focal wave [/+ as a 2D sys- 
tem of homocentric rays, which common point (focus) 
is at the origin (x,y) = (0,0). These rays are re- 
fracted at the plane surface y = —L/2 becoming into 
the field F+, for which the focus is paraxially imaged at 
y = (y/ e c /e — l)L/2 w Re(5)L/4. In our numerical exam- 
ple, the paraxial image would be placed at y = 0.25 jam. 
Also it would exhibit strong spherical aberration. Ap- 
plying a similar picture to the counter-propagating wave 
field U~ we are able to recognize the optical ray tracing 
shown in Fig. [5] for the region near the z axis. In the plot 
we represent for simplicity the projection of ray traces 
onto the xy plane inside the lamellar dielectric structure. 
Our previous discussion supports the high ray concen- 
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Fig. 6. Resultant field e x {x, y) for different Bessel indices 
n in the vicinities of a thin slab with L = 0.5 jam and 6 = 
— 10. Subfigures in the upper row account for amplitudes; 
associated complex arguments are depicted below. 



tration that is produced around the dual optical axes 
(x, y) = (0, ±0.25) jam. We may conclude that Bessel ex- 
citations in the presence of the thin dielectric layer yield 
a pair of localized wave fields whose proximity leads to 
a linear interaction still conserving the diffraction-free 
characteristic of the originating beam. 

In order to analyze surface effects we reduce the slab 
width L approaching a wavelength and, additionally, the 
wave impedance mismatching is increased. In the limit 
\5\ — > oo we obtain = and = — ex.p(—ik y L)U^ 
from Eqs. ([T2j) and ([13]). In this case, the core layer 
behaves like a infinite tunneling barrier prohibiting the 
transit of the wave through itself. Thus each reflected 
field, 

E±(x,y) = -U2(x,-y±L), (15) 

reproduces a mirror-symmetric version of its correspond- 
ing incident wave field in the cladding medium. Also the 
phase factor of argument — k y L in the spatial spectrum 
leads to a displacement L along the y axis of the original 
focus. That is, the geometric focus of the wave field U ± is 
stigmatically imaged by means of R T at (x, y) = (0, =fL), 
respectively. 

Fig. [6] shows the complex amplitude and phase of the 
wave function e x (x,y) given in Eq. (Q3j) for material pa- 
rameters L = 0.5 jam and 5 = — 10. A bifocal 2D pattern 
is clear in the amplitude figures following a previous dis- 
cussion. Interestingly, the size of the twin spots is signif- 
icantly lower in the y direction than that of the source 
fields and it is close to the width of the Jo pat- 
tern. However, this superresolving effect is nearly lost 
if L increases sufficiently. We may conclude that the in- 
terference of propagating and counter-propagating fields, 
whose spectra approaches in y < —L/2 and U~ in 
y > L/2 in the limit \S\ — > oo, yields a couple of intense 
superresolving spots near the surfaces of the thin core 
film. This effect is independent of the Bessel index n, 



n = n = 1 n = -2 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



y (y.m) y (\im) y (\lm) 

Fig. 7. Electric field e x (x,y) given in Eq. (jT4j) for L — 
0.5 jam and 5 = 1. For a given Bessel index n the field 
amplitude is represented above and its complex argument 
below. 

however it drives a lateral shift along the x axis. Finally, 
one can see in Fig. [6] a set of phase singularities, which 
charge is independent of the index n from the Besselian 
source. 

As seen in the amplitude distribution of Figs. [4] and El 
the Bessel index n induces a lateral shift of the twins in 
the presence of the thin film. However, this is not a uni- 
versal rule as shown in the following example. The trans- 
verse profile of the electric field e x (x,y) for L = 0.5 jam 
and 6 = 1 is depicted in Fig. In this case, the wave 
clearly penetrates into the core slab leading to a net 
transfer of energy from one side of the cladding to the 
other side. Then the transmitted field is coupled with 
the counter-propagating incident wave in order to regen- 
erate the Besselian pattern, which however shows a sig- 
nificant distortion and blurring. Inside the core film, the 
transmitted fields are superposed with scattered waves 
originated from multiple reflections on the surfaces. The 
multiple-beam interference leads to a standing wave in 
the xy plane with resonant peaks and nodes (phase sin- 
gularities) mainly along the y axis. For instance, three 
peaks form the shape of the diffraction- free beam if n = 0, 
as shown in Fig. [71(a). The profiles exhibit more complex 
features for higher n orders. In our numerical example, 
however, it gives rise to a pair of conspicuous spots of 
irregular contour if n ^ 0, seen in Figs. [3(c) and[3(e). 

5. Concluding remarks 

In this paper we have reported the formation of highly- 
localized nondiffracting beams in the presence of a core 
microfilm by means of a Besselian source. Differential 
transverse patterns come out consisting mainly of a cou- 
ple of non-centrosymmetric peaks enclosed by numerous 
phase singularities. In general the index n of the Bessel 
function fixes a shift, lying in a plane parallel to the in- 
terfaces, of the center of mass of each spot of light. When 
the slab width is of the order of the wavelength, in prin- 
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ciple, the neighboring twins may be indefinitely propa- 
gated brushing the plane surfaces. Furthermore, a high 
concentration of electromagnetic energy density may be 
additionally found inside the monolayer if the relative 
increment of the dielectric constant is moderate or low. 

Practical implementation of these ideas requires the 
same well-known procedures to produce a Bessel beam. 
For instance, a field spectrum like that of Eq. Q may 
be experimentally attained using an opaque screen, with 
a centered extremely- thin transparent annulus, placed 
at the front focal plane of a perfect lens 26 . Highly- 
efficient approaches may be found using axicon a 31 i 32 , 
Fabry-Perot interferometers 33 , leaky scree mj 34 i 35 i 36 , and 
diffractive optical elements ^ 37 ! 38 . Using such devices as 
external sources in our system would excite the required 
diffraction-free wavefields in the layered medium. Let 
us point out that the experimental conditions may im- 
pose dramatic limitations. For instance, the generation 
of wave packets with mean propagation constant (3 and 
on-axis spectral width A/3 leads to conservation of the 
beam intensity along a finite distance, Az « A/ A/? 26 . 

The theoretical analysis carried out in previous sec- 
tions allows us to tackle more general problems in a plain 
way. If the thin film is deposited on a substrate like glass, 
the ambient medium and the substrate constituting the 
cladding have different dielectric constants. The multi- 
layered optical system is then described by incorporat- 
ing two new wave fields for the substrate region, as 
those given in Eq. ([TO)) for the core layer. In this case 
we would use k 3jV = ^/k 2 — (3 2 — k\ , where k s = uj^Jel/c 
and e s is the dielectric constant of the substrate. Fi- 
nally the spectra F^, ano ^ would be obtained 
from by applying appropriate boundary conditions. 
Also, in practical implementations we might consider a 
7r-illuminating configuration where counter-propagating 
fields do not impinge directly on the substrate. Here we 
set U~ = and wave localization is accomplished solely 
by the source field J7+ travelling in the ambient medium. 
More sophisticated setups might be examined straight- 
forwardly. 

The properties and principles of Bessel-like beams in 
bulk media have offered multiple applications, includ- 
ing precision atomic micromanipulatio n 39 ! 40 ! 41 , laser mi- 
crofabrication of transparent materials 42 , electron laser 
acceleration 43 , and highly-efficient generation of nonlin- 
ear effect o 44 i 45 . The formation of similar light structures 
near surfaces and interfaces leaves an open door to poten- 
tial uses in non-destructive characterization of material 
thin films^ 2 ". Provided grazing incidence onto the core 
slab, such a procedure would increase wave penetration 
and enhance light-matter interaction. 
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List of Figures 

1 Normalized amplitude (upper row) and com- 
plex argument (lower row) of the wave 
function U n (x,y) for n = 9.07 /im" 1 and 
indices (a-b) n = 0, (c-d) n = 1, and (e-f) 
n = -2 

2 (a) Schematic geometry of the planar- 
layer-based medium. In (b) we represent 
the wave fields in regions established by 
interfaces at y = ±L/2. Perfect matching 
e c = e leads to wave fields U ± = F ± = R^. 

3 Wave function U+(x,y) associated with 
Bessel beams of Fig. HJ Similarly, normal- 
ized amplitude is represented in the upper 
row and complex argument is plotted in 
the lower row 

4 Transverse distribution of Bessel-excited 
wave field e x (x,y) inside a thin slab of 
width L — 100 jam made of a material 
with relative dielectric constant 6 — 10 -2 . 
Again \e x \ is plotted in the upper row and 
its complex argument is shown in the lower 
row 
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Ray tracing representation of the Bessel- 
like beams shown in Fig. [H Dual parallel 
axes at y = ±0.25 fim (and x = 0) set 
the confinement region of the nondifiract- 
ing beam. Spherical aberration drawing 
the envelope of the refracted rays gives rise 
to a bilateral-symmetric caustic 



Resultant field e x (x, y) for different Bessel 
indices n in the vicinities of a thin slab 
with L — 0.5 firn and S = —10. Subfig- 
ures in the upper row account for ampli- 
tudes; associated complex arguments are 
depicted below 
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Electric field e x (x, y) given in Eq. ([T4|) for 
L = 0.5 jim and 6=1. For a given Bessel 
index n the field amplitude is represented 
above and its complex argument below. . 



